Energy Dependence of the Cronin Effect from Non-Linear QCD Evolution* 
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The non-linear evolution of dense partonic systems has been suggested as one of the novel physics 
mechanisms relevant to the dynamics of hadron-nucleus and nucleus-nucleus collisions at collider 
energies. Here we study to what extent the description of Cronin enhancement in the framework of 
this non-linear evolution is consistent with the recent observation in -y/s = 200 GeV d-Au collisions 
at the Relativistic Heavy Ion Collider. We solve the Balitsky-Kovchegov (BK) evolution equation 
numerically for several initial conditions encoding Cronin enhancement. We find that the properly 
normalized nuclear gluon distribution is suppressed at all momenta relative to that of a single 
nucleon. Calculating the resulting spectrum of produced gluons in p-A and A-A collisions, we 
establish that the nonlinear QCD evolution is unable to generate a Cronin type enhancement, and 
that it quickly erases any such enhancement which may be present at lower energies. 
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The observation that the ratio of particle yields in p- 
A and A-A, scaled by the number of collisions, exceeds 
unity in an intermediate transverse momentum range of 
a few GeV, is commonly referred to as the Cronin ef- 
fect. This was first seen at lower fixed target energies Q 
and was recently confirmed in ^/s = 200 GeV d-Au col- 
lisions at RHIC 2]. The current interest focuses mainly 
on comparing this Cronin enhancement in d-Au to the 
relative suppression of produced hadrons in Au-Au colli- 
sions at the same center of mass energy and in the same 
transverse momentum range 0- The opposite trend of 
the two effects and their centrality dependence suggests 
that d-Au data may serve as an efficient benchmark mea- 
surement to distinguish between the two different physi- 
cal mechanisms suggested for the relative suppression of 
hadron spectra in Au-Au collisions: initial state parton 
saturation^ and final state jet quenching 

The physics of dense partonic systems and their non- 
linear perturbative evolution to higher energy has mo- 
tivated several attempts at understanding bulk proper- 
ties of ultra-relativistic heavy ion collisions such as the 
multiplicity, rapidity distribution and centrality depen- 
dence of particle production [HQ. In particular it has 
been suggested that saturation effects can account for the 
suppression of the high-p^ hadronic spectra in Au-Au 
collisions at RHIC. On the other hand, it is known that 
saturation models based on multiple scattering (the so 
called Glauber-Mueller Q or McLerran-VenugopalanH 
models) exhibit Cronin enhancement in p-A [ij 0, El 
and A-A 0, . In these models, the saturation of low 
Pt gluons is the result of a redistribution of gluons in 
transverse phase space Q] which does not change 
the total number of gluons, thus resulting in a compen- 
sating enhancement at momenta just above the satura- 
tion momentum Q s . What is not fully understood is i) 



whether such Cronin enhancement encoded in the initial 
condition of a nuclear wave function persists in the non- 
linear perturbative QCD evolution to higher energy and 
ii) whether such Cronin enhancement can be generated 
by the non-linear evolution itself. This paper goes beyond 
earlier discussions [ToL HH by providing the first com- 
plete (numerical) answer to these questions. We do not 
address other approaches to Cronin enhancement |l5| . 

We start from the Balitsky-Kovchegov (BK) evolution 
equation 0, 0], which describes the evolution of the 
forward scattering amplitude N(r, y) of a QCD dipole of 
transverse size |r| with rapidity Y and y = (a s N c /ir) Y , 



dN(\r\,y) 



1 

2^ 



dV (r - z) ' z 



(1) 



dy 2ir J (r — z) 2 z 2 

x [N(\r - z|) + N(\z\) - N(\r\) - N(\r - z\)N(\z\)} . 

The unintegrated gluon distribution is related to the in- 
clusive gluon distribution 4>{k) oc d ^ x ^ x Ji h ^ and is given 
in terms of the dipole amplitude 



d 2 r 
2vrr 2 



exp{i r • k} N(r) . 



(2) 



In the following, we also use the modified gluon distribu- 
tion 



h(k) = k 2 V 2 k <p(k) 



(3) 



*We dedicate this work to the memory of Ian Kogan. 



The two definitions coincide for the leading order per- 
turbative distribution <j>(k) oc -p-, but are different in 
general, and especially at low momenta. 

Using the second order Runge-Kutta algorithm we 
solve the BK equation (TjQ) numerically with 8000 equally 
spaced intervals in lnfc-space between —15 and 35 and 
a step Ay = 0.0025. The accuracy of this algorithm is 
better than 2 % in the entire range of k discussed below. 
We evolve two initial conditions given by the McLerran- 
Venugopalan $ (MV) and Golec-Biernat-Wiisthoff 
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(GBW) model respectively: 
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Q> 2 



Nmv = 1 



exp 



In 
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N, 



Qs 

GBW 



1 — exp 



(4) 
(5) 



where Aq C d = 0.2 GeV. For momenta k > 0(1 GeV), 
the sensitivity on the infrared cut off e is negligible. The 
amplitudes Nmv and Nqbw are similar for momenta of 
order Q s , but differ strongly in their high k behavior; 
4>GBw(k) decays exponentially while <Pmv has a power- 
like tail 1/k 2 . 

Fig^Jshows the evolution of h(k, y) and (f>(k, y) for dif- 
ferent initial conditions. The solutions for h(k, y) quickly 
approach a universal soliton-like shape and do not change 
further except uniformly moving in k on the logarithmic 
plot. The position of the maximum is the evolved value 
of the saturation momentum Q s (y). The solutions for 
different initial conditions and different rapidities scale 
as a function of the scaling variable p = k/Q s (y). The 
shape of the initial condition affects only the value of the 
saturation momentum Q s (y), but not the shape of the 
evolved function h(p,y). The y-dependence of h(p,y) is 
very weak: the function evolves fast towards a scaling 
form h{p). As the rapidity changes between y = 4 and 
y = 10, the ratio h(p, yi)/h(p, 2/2) varies by at most 40% 
over three orders of magnitude of the scaling variable 
p. Similar behavior was found for cj) (results not shown) . 
This is consistent with previous numerical works |20l l2l| . 

To get a quantitative idea of the behavior of the scal- 
ing solution, we fitted the numerical solution of 4>(p) 
to two analytical expressions: si(p) = ap 2 ' 1_A - ) and 
s 2{p) = a\n(bp)p 2 ( 1 ~^ for p > 5. The functional form 
si with A = 0.37 and lnQ s oc y describes the sc aling 
behavior of solutions of the linear BFKL equation [22|. 
It was argued in Ref. [23j that S2 with the same value 
of A and hiQ s oc -j-^j-y — 2(1- A) ^ n 2/ accounts for the 
effects of nonlinearities in Eq. (JJJ. We find that s\ does 
not give an acceptable fit to <fi(p) m any extended range 
of p. For values of p between 1 and 10 3 the value of 
A varies between 0.39 and 0.46. This is in contrast to 
the BFKL equation, where we find numerically that s± 
with A = 0.37 does indeed approximate the solution over 
several orders of magnitude with very good accuracy (re- 
sults not shown). On the other hand, for 5 < p < 1000, 
S2 gives a good fit with A = 0.32. If, following [23| we fix 
A = 0.37, the fit is still good. 

To study the effect of the evolution on the Cronin en- 
hancement, we consider two initial conditions, AnLy(r) 
and N% v {r) with q 2 = 0.1 GeV 2 and Q 2 = 2 GeV 2 . 
Since q ~ ^qcd and Q is of order of the estimated 
saturation momentum for a gold nucleus @], this choice 
mimics the gluon distributions of a proton and of a nu- 
cleus respectively. At large transverse momenta the ratio 
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FIG. 1: Solutions of the BK equation. Upper-left: h(k) 
evolved (left to right) from y = to 5 and 10 for different 
initial conditions: GBW with Q 2 a = 0.36 GeV 2 (solid lines), 
MV with Q 2 a =4 GeV 2 (dashed lines) and MV with Q 2 = 100 
GeV 2 (dotted lines). Upper-right: The same as upper left for 
4>{k). Lower-left: the scaled function h(p) versus p — k/Q s 
for y = 4, 6, 8, 10 and the same initial conditions and conven- 
tions (lines cannot be distinguished). Lower-right: Ratio of 

Ky,p)IKv;p = 1) over Kv = w,p)/Hy = 10 >p = !) for 

y — 4 (solid line), 6 (dashed line), 8 (dotted line) and 10 
(dashed-dotted line), and initial condition MV with Q 2 S = 4 
GeV 2 . 



of the corresponding Fourier transforms is given by the 
ratio of the saturation momenta, 



fe°(fc,y = Q) 

hi(k,y = 0) 



91 = 4 1/3 

2 
q- 



(6) 



This relation also holds for <j>. As discussed in [ToL Hl| , 
these initial conditions exhibit Cronin enhancement, 



namely 



h Q (k,y=0) 



> 1 for k ~ Q. We solve the BK 



A 1 / 3 hi(k,y=0) 

equation with these two initial conditions and construct 
the ratio R(k,y) = h c ^(k,y)/A 1 ^ 3 h q (k,y) and the corre- 
sponding ratio for 0, see Fig. |2 The initial Cronin en- 
hancement at rapidity y = is seen to be wiped out very 
quickly by the evolution. Within less than half a unit 
of rapidity y the ratios show uniform suppression for all 
values of transverse momentum. The observed behaviour 
persists if different amounts of Cronin enhancement arc 
included in the initial condition. 

As seen in the lower panel of Fig. |3 the Cronin en- 
hancement also disappears rapidly with rapidity when 
gluon distributions are evolved according to the linear 
BFKL equation. Qualitative differences between the 
BFKL and BK dynamics are only visible at momenta 
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FIG. 2: Ratio of distributions (j> and h in nucleus and proton, 
normalized to 1 at k — > oo. Upper plots: BK evolution, with 
MV as initial condition with Q 2 S = 0.1 GeV 2 for p and 2 GeV 2 
for A. Lines from top to bottom correspond to y = 0, 0.05, 
0.1, 0.2, 0.4, 0.6, 1, 1.4 and 2. Lower plots: BFKL evolution, 
with MV as initial condition with Q\ — 4 GeV 2 for p and 100 
GeV 2 for A. Lines from top to bottom correspond to y = 0, 
1 and 4. 



k < Q s , where saturation effects are important. For 
larger momenta k, the ratios are very similar for linear 
and non-linear QCD evolution. We thus conclude that 
the wiping out of the initial enhancement is primarily 
driven by the linear BFKL dynamics which is contained 
in the BK equation as well. 

For the evolved gluon distributions determined above, 
we have calculated the yield of produced gluons in p-A 
and A-A collisions at central rapidity according to the 
factorized expressions 24] 

oci/ d 2 kh"(y,k)h^(y,p-k), (7) 
aya z pa z b p z J 

a!N AA oc —-TT- [ d 2 khQ(y,k)hQ(y,p-k). (8) 



dyd 2 pd 2 b p 2 

From these spectra, we calculate the p- and y-dependent 
ratios 



RpA 



dN pA 
dyd 2 p d 2 b 

dyd 2 p d 2 b 



dN A j 



Raa = 



dyd 2 p d 2 b 
dyd 2 p d 2 b 



As seen in Fig. |3 the non-linear BK evolution quickly 
wipes out any initial Cronin enhancement not only on 
the level of single parton distribution functions but also 
on the level of particle spectra. We have checked that 
this behaviour is generic by evolving different initial con- 
ditions corresponding to different initial amounts of en- 



hancement. We note that in calculations of the gluon pro- 
duction in p-A in the eikonal approximation |l3|, |2fj|, |26( , 
the gluon distribution h rather than tj> enters the right 
hand side of Q, but no similar statement exists of 
nucleus-nucleus collisions. We have checked that the re- 
sults using 4> are very close to those shown in Fig.|3J More 
generally, the expressions J7J) and JSJl are based on rather 
strong approximations discussed in Rcf. 10] . However, 
our conclusion about the disappearance of Cronin en- 
hancement during QCD evolution is likely to persist in 
more refined ways of calculating particle spectra, since 
it is rooted directly in the rapidity dependence of gluon 
distributions. 
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FIG. 3: Ratios R p a and Raa of gluon yields in p-A (upper 
plot) and A-A (lower plot) for BK evolution, with MV as 
initial condition with Q 2 =0.1 GeV 2 for p and 2 GeV 2 for A. 
Lines from top to bottom correspond to y = 0, 0.05, 0.1, 0.2, 
0.4, 0.6, 1, 1.4 and 2. 

We now comment on a recent formal argument [ill ] 
which - in contrast to our numerical findings - suggests 
that Cronin enhancement survives the non-linear evolu- 
tion. It is based on the observation that at very short dis- 
tances r — > 0, the dipole amplitude N(r) is not affected 
by evolution. Thus, the integral of the gluon distribution 
function <j) over the transverse momentum is expected to 
be rapidity independent, 



d 2 k(j>{k) 



1 



N(r)\ r 



(9) 



One thus obtains the sum rule 
d 2 k<f> A (k,y) = A 1 '* 



d 2 k(j) p (k,y) (10) 



valid for any rapidity, since it is satisfied by the initial 
condition <Pmv- Since the nonlinear evolution leads to 



4 



the depletion of the gluon distribution 4>a (k) relative to 
A 1 ^ 3 (j) p (k) at low momenta, it must follow that in some 
range of momenta this effect is compensated by enhance- 
ment of 4>A- Similar arguments can be made about the 
distribution h and also about the gluon yield in p-A and 
A- A collisions. 

However, this argument breaks down since the quan- 
tity defined in Eq. © is infinite. As such Eq. (fTT)f> relates 
only the divergent parts of the integrals which are dom- 
inated by ultraviolet, and carries no information about 
possible behaviour at finite momentum. To be more spe- 
cific, we use the scaling property 4>(k,y) = (f>(k/Q s (y)) 
of the solution of the BK equation established above. 
It is known that the ratio of the saturation momenta 
for any two solutions is preserved by the BK evolution 
[20I I21I |22| . For our two solutions representing a nu- 
cleus and a nucleoli, this implies qp = A 1 / 6 . We now 
rewrite the sum rule H1(J|) by regulating the divergent in- 
tegrals with a large but finite UV cutoff aQf(y), 

f a 2 (Qf) 2 f a 2 

/ d 2 kMk,y) = (Qf) 2 d 2 P <p(p) (11) 

Jo Jo 

ra pa (QJ) 

= A^iQP) 2 / d 2 p^{p) = A 1 / 3 / d 2 kMk, y) ■ 
Jo Jo 

In the formal limit a — > 00, we recover Eq. (llUfl . How- 
ever, since Qf ^> Q p s , the regularized sum rule 111|) is 
easily satisfied even if the nuclear distribution is sup- 
pressed relative to that of a single nucleon uniformly at 
all momenta. Thus, the sum rule H10|) carries no infor- 
mation about either presence or absence of the Cronin 
enhancement. 

In summary, we have found that the nonlinear QCD 
evolution to high energy is very efficient in erasing any 
Cronin type enhancement which may be present in the 
initial conditions. For realistic initial conditions this dis- 
appearance occurs within half a unit of rapidity. We 
note that in our units the evolution from 130 GeV to 
200 GeV corresponds to Sy ~ 0.1 for a s = 0.2, and 
thus is not sufficient to completely eliminate an initial 
enhancement at central rapidity. For forward rapidity, 
6y is greater. The evolution to the LHC energy corre- 
sponds to Sy ~ 1. Thus the BK evolution suggests the 
reduction of the Cronin effect in d-Au for forward ra- 
pidities at RHIC and it predicts the disappearance of the 
Cronin effect for p-A collisions at LHC. 
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